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Abstract 



We outline a method for calculating the self force (the "radiation reaction 
force") acting on a scalar particle in a strong field orbit in a Schwarzschild 
spacetime. In this method, the contribution to the self force associated with 
each multipole mode of the particle's field is calculated separately, and the 
sum over modes is then evaluated, subject to a certain regularization pro- 
cedure. We present some explicit results concerning the implementation of 
the calculation scheme for a static particle, and also for a uniform circular 
motion, on the Schwarzschild background. 
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This short manuscript deals with the calculation of the self force acting on a small charged 
object (a "particle") moving in a Schwarzschild spacetime. We shall briefly outline here a 
method, recently proposed by one of us (A.O.) and present some results. This method 
involves the calculation of the contribution to the self-force due to each individual field's 
multipole, and the summation over these contributions, subject to a certain regularization 
procedure (outlined below). For simplicity we shall consider here a scalar charge. We 
expect, however, that the generalization to the electromagnetic self force will be almost 
straightforward. The generalization to the Kerr case and to the gravitational self force is 
likely to be possible as well. A closely related approach to the regularization of the self 
force, also based on the multipole expansion, is currently being studied by Lousto 0. For 
an overview of further related works, the reader may refer, e.g., to Ref. P]. 

The main challenge in calculating the self force is to deal with the various infinite quan- 
tities associated with the divergence of the field at the particle's location. Previously it 
has been proposed |^ that radiation reaction effects on the orbit may be calculated by 
evaluating the contribution of each Fourier-multipole mode imw of the field to the radiative 
evolution, and then summing over all modes. This approach has two advantages: First, each 
individual mode of the field turns out to be continuous (and the corresponding self force to 
be bounded) even at the particle's location. Secondly, calculating each Imw mode of the 
field becomes a relatively simple task, as it only requires the solution of an ODE. In Ref. 
@] this sum-over-modes approach has been proposed for the calculation of the adiabatic, 
orbit-integrated, evolution rate of the three constants of motion in Kerr: The Energy E, the 
azimuthal angular momentum L, and the Carter constant Q. One might have hoped that 
the same method could also be used for calculating the momentary self force. Unfortunately, 
at least for the self force, this naive procedure turns out, in general, inapplicable: Although 
each mode yields a finite contribution, the sum over all modes is found to diverge. This is 
the situation even in the simple case of a static particle in a fiat space: For such a particle, 
located at a distance tq from the origin of coordinates (with respect to which the spherical 
harmonics are defined), the contribution of each multipole ^ to the radial component of the 

self force is the same g^/(2rQ), where q is the particle's scalar charge. Obviously, the 

sum over £ diverges.^ The regularization scheme outlined here is aimed to overcome this 
type of divergence (though in a nontrivial black- hole spacetime). 

Practically, in this approach the calculation of the self force is composed of two separate 
parts: 

(i) The numerical part: solving the appropriate ODE for each mode imw of the field, and 
evaluating each mode's contribution to the self forceQ — a quantity we denote by F^"^'^. 
(Alternatively, one may numerically solve the 1+1 PDE in the time domain, for each £ and 
m.) 



Note, however, that the mode sum for the adiabatic, orbit-integrated, evolution rate of E and L 
does converge [Q. It is not clear yet whether the corresponding mode sum for Q converges or not. 

^We use here = qcj)^^ as the basic expression for the force applied on a charged particle by the 
scalar field. Throughout this manuscript we shall often omit the vectorial index of the force and 
related quantitied for notational brevity. 
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(ii) The analytic part: applying a certain regularization procedure to the mode sum. This 
procedure involves the analytic calculation of certain regularization parameters, as we de- 
scribe below. 

This manuscript deals with the second part, namely, the analytic regularization scheme. 
The first part — the numerical determination of F^"^'^ — is currently being implemented by 
Burko for several physical scenarios 

The analytic regularization scheme has already been fully implemented by one of us 
(A.O.) for circular orbits in the Schwarzschild geometry [jlO[. Barack and Ori are currently 



developing the scheme for noncircular orbits (focusing on radial orbits at the moment). Since 
our regularization method has already been implemented in practical calculations (for static 
P and circular P orbits in the Schwarzschild geometry), we found it useful to provide a 
short account of the regularization scheme, and to briefly describe the main results already 
obtained. A more detailed account, including a systematic development of the method and 
a detailed calculation of the various parameters involved, will be given in Ref. 

As was shown by Quinn and Wald (see also [|TI| , p!2| ), the physical self force is (in 



vacuum) the sum of two parts: (i) A local, Abraham-Lorentz-Dirac (ALD)-like term, and 
(ii) a "tail" term associated with the tail part of the Green's function. The local 

term is trivial to calculate. We shall therefore focus here on the tail term. This term may 
be expressed as 

i7^{tail) ^ F (1) 

where F^ denotes the contribution to the force (evaluated at r = 0) from the part r < — e 
of the particle's worldline. Decomposing F^ into £-modes, one finds 

^(taii) ^ j.^^ J2 = lim^ J2 - ^Pe) ■ (2) 

Here, F/, and F^ denote the force associated with the £-multipole of the field, sourced, 
respectively, by the interval r < — e, the interval r > — e, and the entire world-line. The 
quantity F^ may be identified with the sum over m, oj of the contributions from all modes 
£, m, uo for a given ^. (Recall that in calculating a field's mode m, uj one takes the source 
term to be the entire worldline. Since we are using the retarded Green's function, the 
part r > does not contribute. However, the interval from — e to 0+ does contribute. 
Essentially, it is this part which is responsible to the instantaneous, divergent, piece of the 
Green's function, which should be removed from the expression for F'-*^''^) 

A clarification is required here concerning the meaning of the last equality in Eq. (^: 
Let To denote the value of r at the self force evaluation point. Then, the quantity F^ is well 
defined at r = tq. The situation with F^ and 5F^ is more delicate, however. Each of these 
quantities has a well defined value at the limit r ^ r^" , and a well defined value at the limit 
r — Tq". Generically, for the r-component (and in some cases also for other components) 
these two one-sided values are not the same. Eq. (0) should thus be viewed as an equation 
for either the limit r ^ oi the relevant quantities (i.e. F^ and 5-F/), or the limit r ^ 
of these quantities. Obviously, this equation is also valid for the averaged force, i.e. the 
average of these two one-sided values. In what follows we shall always consider the averaged 
force. [Of course, the final result of the calculation, F'^^'^^^^ (which has a well defined value 
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at the evaluation point), should be the same regardless of whether it is derived from its 
one-sided limit r ^ rQ , or from r — >■ r^, or from their average.] 

Next, we seek an e-independent function h^, such that the series J^ei^"^ " h^) converges. 
Once such a function is found, then Eq. (Q) becomes 

jT^Ctaii) = Y^[f'- h') - lim^ {^^e - h') . (3) 

In principle, can be found by investigating the asymptotic behavior of as i ^ oo. It is 
also possible, however, to derive from the large-£ asymptotic behavior of 6F^ (the latter 
and F^ must have the same large-£ asymptotic behavior, because their difference yields a 
convergent sum over i). In addition to h^, the investigation of 6F^ should also produce the 
parameter 




which is required for the calculation of F*-*^''^ in Eq. (^) . 

From Eqs. @ and (|) it is obvious that we only need the asymptotic behavior of SF^ 
in the immediate neighborhood of e = 0. It is therefore possible to analyze 6F^ using 
local analytic methods. In particular, we can apply a perturbation analysis to the £-mode 
field equation (in the time domain). To that end, we express the £-mode effective potential 
V^'{r) as a small perturbation 6V^{r) over the value of V^{r) at the evaluation point, Vq = 
V^{r = ro). This perturbative approach is used to analyze = G^fx", (r)] (the £- 
mode Green's function), where x°' is the evaluation point and Xpi^r) describes the particle's 
worldline. Expressing as a function of r and z = ri, the perturbation analysis provides an 
expression for G^ as a power series in r (with 2;-dependent coefficients). Only the first three 
terms, i.e. the terms up to order r^, are required for the calculation of the self force (recall 
that eventually we take the limit e — > 0; as it turns out, in this limit the contributions from 
all higher-order terms vanish). The perturbation analysis yields explicit analytic expressions 
for the required three expansion coefficients (as functions of z), which we present in Ref. 0. 

We next construct 6F^ from G^ (essentially by integrating the latter's gradient from 
r = — e to r = 0). One then obtains the large-^ asymptotic form 6F^ = ai+b+ci~^ +0{i~'^), 
in which the parameters a, b, and c are independent of i and e. (These parameters do depend 
on the evaluation point, and also on the orbit - through the values of the particle's four- 
velocity and its first few proper-time derivatives there.) The regularization function thus 
takes the form 

= ai + b + cr\ (5) 
Consequently, the tail part of the self force is found to be 

^(taii) = Y^(^pi 1) - d. (6) 

e 

The implementation of this regularization scheme thus amounts to analytically calculating 
the four regularization parameters a, b, c, and d, using the above mentioned perturbation 
expansion for the Green's function. (Note that there exist four such parameters for each 
vectorial component of the force.) 
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We shall now present the main results already obtained for the values of the above 
regularization parameters in various cases. These results concern the r-component of the 
self force, and are valid for either geodesic or non-geodesic orbits. It can be shown that the 



parameter a vanishes for all orbits in Schwarzschild [1^, and preliminary results indicate 
that the same holds for c. (Recall that we consider here the averaged force: the "one-sided" 
values of a are, in general, nonzero 0.) On the other hand, the parameters h and d are 
generically nonvanishing. In what follows we give the values of these parameters for some 
physical scenarios. Q 

For static or circular orbits the parameter d (as well as a and c) is found to vanish [|TU|,^: 

ttr = Cr = dr = Q for static or circular orbits, (7) 

where the sub-index r indicates that these parameters are associated with the covariant 
r-component of the self force. Hence, in these cases one finds0 

^jtaii) ^ ^ (^pt. _ {^iaiic or circular orbits). (8) 

For a static particle, the parameter hj. is given by [§] 

f ro- M \ 



— j (static particle), (9) 



where M is the black hole's mass. For a circular orbit, br takes the form [|T0| 

1 / ^ ro-3M 
'2r§ I " ro-2M 



= -TT^ 2 4 - _ h ) (circular orbit), (10) 



where 



4 = F(l/2,l/2;l;t;2) , J, = F(l/2, 3/2; 1; t;^), (11) 



F denotes the hypergeometric function, and v is the tangential velocity with respect to a 
static observer: v = ^ (for an equatorial orbit). 

Recently, Burko numerically calculated F^"^'^ for static and circular [|| orbits, and 
used the above regularization scheme to calculate the self force. The above expressions (PiP^Of) 



for br (as well as the vanishing of and Cr) are in excellent agreement with the large-£ limit 
of his numerically-deduced F^ (obtained by summing F^"^'^ over m and u). 

For radial motion in flat space (i.e., motion with fixed angular coordinates 9 and (p), the 
parameter br is found to be ^ 

br = 7^ (l — r"^ + rgr ) (radial motion in Minkowski), (12) 

2ro ^ ^ 



^We use metric signature ( — h ++) throughout. 

^For static or circular orbits, the r-component of the local, ALD-Uke, term of the total self force 
vanishes. Thus, Eq. (^) describes, if fact, the (r-component of the) total self force. 
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where a dot denotes differentiation with respect to proper time. In this case, the parameter 
d does not vanish, in general, and one finds 

dr = -q'^ (^dr — UrO^^ (radial motion in Minkowski), (13) 

where Mq, and Oq, are the four-velocity and four-acceleration, respectively, and = a'^aa- 
Remarkably, dr thus equals the r-component of the ALD-like local term (given in Ref. 0). 
As a consequence, this contribution from the tail term will cancel out the r-component of 
the local ALD-like term in the expression for the total self force. 

Generalization of the above results, Eqs. (|T2]) and (p!B|), to general radial motion in a 



Schwarzschild spacetime (and also in a more general class of static spherically symmetric 
spacetimes) will be presented in a forthcoming paper 
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